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Abstract. I describe how to obtain a rather good experimental determination of the 
eccentricity of the Earth orbit, as well as the obliquity of the Earth rotation axis, by measuring, 
over the course of a year, the elevation of the Sun as a function of time during a day. With 
a very simple "instrument" consisting of an elementary sundial, first-year students can carry 
out an appealing measurement programme, learn important concepts in experimental physics, 
see concrete applications of kinematics and changes of reference frames, and benefit from a 
hands-on introduction to astronomy. 
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1. Introduction 

One of the cornerstones of introductory courses in classical mechanics is the derivation of 
Kepler's laws. In particular, the derivation of Kepler's first law, stating that the trajectory 
of a planet is an ellipse with the Sun located at one of the foci, is an important application 
of Newton's laws to a multidimensional problem. However, very few students are aware of 
the fact that the eccentricities of the planets of the Solar system are actually quite small, with 
trajectories very close to a circle, which makes Kepler's achievement (based on Tycho Brahe's 
measurements) even more remarkable. 

Here, I describe a simple measurement programme, suitable for first-year university 
students, consisting in measuring the elevation of the Sun as a function of time during a 
day, and in repeating this typically once a week over a full year. By measuring the maximal 
elevation h max of the Sun, and the time t max at which this maximum occurs (i.e. the true 
local noon), students can readily check that these quantities vary a lot over the year. The 
change in /i max is essentially related to the obliquity e of the Earth over the ecliptic, and thus 
allows for quite an accurate determination of e (as well as that of the latitude of observation). 
The change of i max over a year gives an experimental determination of the equation of time 
E{t), i.e. the difference between the mean local noon and the true local noon, and allows 
for a determination of the eccentricity e of the Earth orbit [1]. This is a rewarding result for 
students to realize that with such simple measurements they can obtain good experimental 
values for the above quantities, and that with careful observations one can perform 'science 
without instruments' as did the astronomers of various antique civilizations [2, 3], 

This article is organized as follows. I first describe how to measure in a simple way the 
elevation of the Sun versus time over a day, with an accuracy of about 1°. Then I give the 
results I obtained for h max (t) and E(t) by repeating the measurement about once a week for 
one year, starting in August 2010. I show how one can extract the obliquity e of the Earth's 
axis and the eccentricity e of its orbit by fitting the experimental data with simple, analytic 
expressions. Finally, possible extensions of the work are proposed. Appendix A contains a 
brief reminder on basic notions of spherical astronomy, and should be read first by readers 
not familiar with these notions. In the remaining appendices, the derivation of the analytic 
expressions used for fitting the data is given, so that the article is self-contained. 

2. Measurements 

We are interested in studying the motion of the Earth around the Sun. Using the relativity of 
motion, we can thus simply measure the apparent motion of the Sun on the celestial sphere, 
i.e. the time-dependance of two angles that define the position of the Sun in the sky. 

As we shall see, for our purpose, it is sufficient to measure the elevation of the Sun (also 
called altitude, or height) above the horizon, i.e. the angle h shown in Fig. 1. This can be 
done very simply by measuring the length s of the shadow of a vertical gnomon (i.e. a rod 
with a sharp point) of length I. Then the elevation of the Sun is given by h = arctan(£/s). 

Contrary to the case where one would measure also the azimuthal position of the Sun, 
here, the orientation of the horizontal base does not need to be fixed. One of the advantages 
of using such a simple setup is therefore that one can change the position of the sundial over 
the course of the day, e.g. in order to operate indoor. 
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Figure 1. Schematic view of the elementary sundial. 



2.1. Construction and use of an elementary sundial 

In practice, I used as a gnomon a steel nail protruding from a plywood base of size 20 x 20 cm 2 . 
In order to have the nail as orthogonal to the base as possible, a hole with a diameter slightly 
less than that of the nail was first drilled into the plate using a drill press. The nail I used has a 
length I — 69 mm above the plate. To measure h, one simply installs this elementary sundial 
on a horizontal surface in the sunlight, and measures with a ruler the length of the shadow. 
Two effects limit the accuracy of the measurement: first, due to the final angular diameter of 
the Sun, the shadow is slightly blurred; second, the horizontality of the base when installed on 
the floor of a room, or on a table, is not perfect [4]. In practice, an accuracy of typically one 
degree is easily obtained (this can be estimated quantitatively by repeating the measurement 
several times with the sundial in different positions, in a short interval over which h barely 
varies, and observing the dispersion of the results). 

Concerning the determination of the time t at which h(t) is measured, an accuracy 
of one minute is sufficient for our purpose, and thus a simple wristwatch can be used. 
However it is wise to check that the watch indicates the correct time before starting a series 
of measurements. Nowadays, this can be done very easily using the websites of national time 
agencies [5] that give access to the legal time with an accuracy of one second or better. 

2.2. Measuring the altitude of the Sun over a day 

Figure 2 shows two measurements of h(t), where t is the legal time, performed in Toulouse, 
France (latitude (p = 43° 36' N, longitude A = 1° 27' E) at two different dates. It is very clear 
from the data that the maximal height h mscK of the Sun depends on the date; this is in general 
well known as it is related to the cycle of seasons. 

However, what appears also clearly on Fig. 2 is that the time £ max at which this maximum 
occurs also depends on the date; this however is not widely known by the students, nor even 
by some physicists, probably because the effect is relatively small (a few minutes) though 
perfectly measurable even with our crude setup. 

In order to proceed, we need to extract from hit) the two quantities h max and i max - The 
theoretical expression of h(t) is derived in Appendix B; however we can at this stage keep 
an empiric approach and just fit the data with a simple function. As h(t) is symmetric about 
t = imax (provided one neglects the motion of the Sun with respect to the fixed stars over 
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Figure 2. Measured elevation of the Sun over the course of a day, for two different dates. The 
solid lines are fits by a polynomial (see text). The dashed lines display the values of /i max and 

^max ■ 



a few hours, which is reasonable given the accuracy of our measurements), I chose to fit the 
data with the following polynomial: 

3 



where the five adjustable parameters are h max , t max , and the coefficients h 2 ,4,6- I chose to go 
up to sixth order in order to get a nicer fit at small elevations (in the mornings and evenings) 
but if the data is taken only for a few hours around i max (±3 to 4 hours around i max are 
enough to determine the quantities of interest) one can use only a fourth-order polynomial 
without affecting the results. Such fits are shown as solid lines in Fig. 2. 

The accuracy in the determination of h max and i max obviously depends on the number 
of measurement points; for the data presented in Fig. 2, they are respectively of about 0.2° 
and 1 min, as data points were collected for several hours before and after t max , at a rate of 
four points per hour typically. When the weather is partly cloudy, one sometimes has to stop 
taking data for a while, and the accuracy in the determination of h maK and f max is thus not as 
good. 

2.3. Annual variation of hmax' determination of Lp and e 

By repeating the above measurements typically once per week for a year, the annual variations 
of /i m ax and t max can be obtained. Figure 3 shows the maximal elevation /i max of the Sun as 
a function of time. One observes a quasi-sinusoidal variation with a period of one year. 

The qualitative explanation for this phenomenon is simple, and is usually part of 
geography courses in elementary education, but I repeat it here for completeness. Due to 
the obliquity e of the Earth rotation axis, the angle between the Earth axis and the line joining 
the center of the Earth to the Sun varies between ir/2 — e (at the summer solstice, when the 
Earth axis leans towards the Sun) and tt/2 + e (at the winter solstice). Correspondingly, the 
angular distance between the Sun and the celestial pole varies with a period of one year, and an 
amplitude of e. This induces a similar variation of h max , around a mean value which depends 
on the observer's latitude. For a quantitative treatment, the reader is referred to Appendix C. 
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Figure 3. Measured maximal elevation of the Sun (points with error bars). We observe a 
sinusoidal variation around 7r/2 — tp, with and amplitude of 2e. The solid line is a fit to the 
simple model discussed in the text. 

2.4. Annual variation oft max : determination of e 

We now turn to a more subtle measurement concerning the variation of i max , which defines 
the true local noon. It is convenient to convert the measured values into a quantity called the 
equation of time, that we shall denote by E, defined as the difference between the mean local 
noon and the true local noon (our measured i max ). The former is obtained from the legal 
noon, given by clocks (corrected if necessary by one hour in summer due to daylight saving 
time) by adding (subtracting) four minutes for each degree of longitude west (east) from 
the reference meridian of the corresponding time zone. For instance, in Toulouse (longitude 
A = 1° 27' E), one needs to subtract i max from 13.00 h in winter time and 14.00 h in summer 
time, and then subtract another 5.8 min to correct for the longitude, to obtain the equation 
of time E. For instance, on October 17th (see Fig. 2), we have i max = 13.66 h, thus giving 
E = 60(14 - 13.66) - 5.8 = 14.6 min. 

Figure 4 gives the results obtained by measuring E(t) over a year. One observes a non- 
trivial behavior, the equation of time varying between a maximum of about 16 min in autumn 
and a minimum of about —15 min in winter, and vanishing at four different dates. 

Physically, the origin of the equation of time lies in the fact that the duration of the true 
solar day, i.e. the time elapsed between two successive transits of the Sun across the observer's 
meridian, is not constant over a year. The solar day would have a constant duration if, along 
the year, the Sun moved on the celestial sphere (i) at constant angular velocity, and (ii) along 
the celestial equator (this defines the so-called mean Sun; the time between two transits of 
the mean Sun defines the mean solar day of 86,400 s). However, these two assumptions are 
both wrong: since the Sun moves along the ecliptic, which is inclined with respect to the 
equator due to the obliquity e of the Earth axis, the motion of its projection on the equator is 
irregular (it coincides with the mean Sun at the time of the equinox, then lags behind the mean 
Sun for a quarter of the year, catches up at the solstice, and then is ahead of the mean Sun 
for another three months). This contribution E\ to the equation of time thus has a 6-month 
period. Moreover, via Kepler's second law of areal velocity (see Appendix D) the angular 
velocity of the apparent motion of the Sun is not constant over the year, due to the fact that the 
Earth orbit is not circular: for instance, when the Earth-Sun distance is smaller (in January), 




Figure 4. (a) Measured equation of time E (points with error bars). The solid line is a fit to the 
simple model of the equation of time given in the appendix, (b) Combining the effect of the 
Earth obliquity (E\(t), dotted line) and of the eccentricity of the Earth orbit (E2(t), dashed 
line) gives the full expression of the equation of time E(t) (solid line). 



the Sun moves faster along the ecliptic. This contribution E^ to E(t) has obviously a one-year 
period. Combining these two contributions explains the temporal variation of the equation of 
time (see Fig. 4b; Appendix D gives the derivation of the analytical expressions of E\ and 
E 2 ). 



3. Exploiting the data 

3.1. Obliquity of the Earth 

It is easy to show (see Appendix C) that to a very good approximation, one has: 

h max = - - ip + e sin f — (t - t ) j , (2) 

where tp is the latitude of the place of observation, e the obliquity of the Earth axis, T 
the duration of the year, and to the date of the vernal equinox. When fitting the data by 
equation (2) with the four previous quantities as adjustable parameters, we obtain 

tp = 43.8 ±0.2°, 
e = 23.5 ±0.1°, 

T = 374 ± 6 d, ( ' 

t = 78 ± 1 d, i.e. March 19 th , 

which is close to the accepted values (respectively, 43.60°, 23.44°, 365.25 d, and March 
20 th .) Note that by repeating the measurements over the course of several years, a much more 
accurate determination of the duration T of the year could be achieved. 



3.2. Eccentricity of the Earth orbit 

We show in Appendix D that a good approximation of the equation of time is given by: 
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where d is the duration of a day (i.e. 1440 minutes), e the eccentricity of the Earth orbit, and 
ti the date of perihelion passage. 

Fitting the data shown in Fig. 4 by equation (4) with e and t\ as adjustable parameters 
(and using the values determined above for e and T), we obtain 



again in relatively good agreement with the values e = 0.0167 and t\ = 3 d found in the 
literature. 

4. Conclusion and outlook 

I have shown that with very modest equipment, one can measure with reasonable accuracy 
some of the orbital elements of the Earth, and in particular its eccentricity, despite its relatively 
small value. The above measurements can be the basis of further activities for students. 
Among them, one can list the following ones, given here under the form of exercises: 

• Use Eq. (B.4) of appendix Appendix B to calculate the length of daytime as a function 
of the latitude along the year and compare it to the one obtained from the ephemerides 
given in calendars. 

• Show that the duration of a solar day is (1 + dE/dt) x 86,400 s. What are its minimal 
and maximal values? 

• Using a sundial with a fixed base, check experimentally that the azimuthal position of 
the Sun when reaches its highest elevation is always the same (i.e., South) throughout 
the year [6]. 

• Still with a fixed-base sundial, plot experimentally the curve traced out over the year 
by the end of the shadow at the mean local noon. This eight-shaped curve is called an 
analema and is sometimes encountered on sundials in order to allow for a computation 
of the legal time from the measured solar time. 
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Appendix A. A quick reminder on spherical astronomy 

This section consists in a minimalist reminder about basic terms and notions of spherical 
astronomy needed for the understanding of the paper. The reader is referred to the first 
sections of Ref. [2] for a similar but more exhaustive reminder. A very detailed and accessible 
introduction to spherical astronomy can be found online in the celestial mechanics lecture 
notes of Ref. [7], or in standard textbooks about spherical astronomy [8]. 

To an Earth-bound observer O, celestial bodies appear to move on a sphere centered on 
himself, the celestial sphere (Fig. Ala). The local vertical points toward the zenith Z; the 
great circle perpendicular to the vertical is the horizon. Over a day, "fixed" stars appear to 
rotate around the north celestial pole P (close to the star Polaris). The great circle going 
through Z and P crosses the horizon in two points defining the South S and the North N; the 



{ 



e = 0.017 ±0.001, 

t\ = 1 ± 5 d, i.e. January 1 st , 



(5) 
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Figure Al. Geometry of the celestial sphere, (a): important points, axes and circles on the 
celestial sphere, (b): definition of the coordinates in the horizon system, (c): definition of the 
coordinates in the equatorial system, (d): definition of the coordinates of the Sun in the 
equatorial system during its motion along the ecliptic. 

two other cardinal points on the horizon (East E and West W) are deduced from S and N by 
a 90° rotation around the vertical axis. 

The position of celestial bodies can be specified in spherical coordinates by two angles, 
once a reference system has been chosen. Two reference systems are particularly useful: the 
horizon system, in which measurements are made in practice, and which is dependent on the 
observer's location on Earth, and the equatorial system, which is defined by the directions of 
"fixed" stars. 

• In the horizon system (Fig. Alb), the local vertical, pointing towards the zenith, is chosen 
as the z-axis. The position of a point on the celestial sphere is specified using the altitude 
(or elevation) h (angle between the radius-vector of the point and the horizon), and 
the azimuth rp, counted along the horizon, starting from the South and counted positive 
towards the West (note that other conventions exist for the choice of the azimuth origin). 

• In the equatorial system (Fig. Ale), OP is chosen as the polar axis. The intersection of 
the plane perpendicular to OP with the celestial sphere defines the celestial equator. The 
position of a point on the celestial sphere is given by the declination S (angle from the 
celestial equator to the body) and the hour angle H (counted along the celestial equator, 
from the south to the equatorial projection of the body). 

The equatorial system is obtained by rotating the horizon system around OE by an angle 
7r/2 — ip, where <p is the (geographical) latitude of the point of observation; for instance, in at 
the North pole (<p — ir/2) the equatorial and horizon systems coincide. 
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Finally, the ecliptic is the great circle along which the apparent annual motion of the Sun 
(traditionally denoted by the astronomical symbol 0) takes place on the celestial sphere. It is 
inclined on the celestial equator by the obliquity e of the Earth axis. 

The angular position of the Sun along the ecliptic is given by the ecliptic longitude A©, 
whose origin is taken at the vernal point V (the point where the ecliptic crosses the celestial 
equator, and where the Sun is located at the time of the spring equinox in the Northern 
hemisphere). The angular distance on the celestial equator between V and the projection 
of the Sun on the equator is the Sun's right ascension a©. 



Appendix B. Expression of h(t) over a day 

From the definitions given above, in the equatorial coordinate system (x',y',z r ), the 
coordinates of the Sun read: 

(x' = cos S cos H \ 
y' = -cos<5 smH , (B.l) 
z' = sin S J 

with 6 the declination and H the hour angle of the Sun. In the horizon system (x, y, z), they 
read 

(x = cos h cos ip \ 
y = — cos h sin tp , (B.2) 
z — sin/i I 

where h is the altitude of the Sun, and tp its azimuth. Since the equatorial system is deduced 
from the horizon system by a rotation of angle ir/2 — ip around the (Oy) = (Oy') axis, the 
(x, y, z) coordinates are obtained by multiplying the (x' , y' , z') ones by the following rotation 
matrix: 

(sin ip — cos ip \ 
1 . (B.3) 

cos (p sin ip I 

From the last component of the relation S = RS' we get 

sin h = cos ip cos 6 cos H + sin <p sin 8 (B.4) 

which gives the elevation of the Sun as a function of time (i.e. the hour angle H) for given 
location and declination of the Sun. 



Appendix C. A simple model for h max (t) 

From equation (B.4) above, we find immediately that the elevation of the Sun becomes 
maximal when cos H = 1 and reaches the value /i max fulfilling 

sinft, max = coscp cos 5 + simp sin S = cos(p — S) (C.l) 

whence 

TT 

h max = — - (p + S. (C.2) 

We now need to express the time dependance of the declination of the Sun. Using Fig. Al, 
one can show (see Eq. (D.5) below) that sin(5 Q = sine sin A©, which can be simplified to 
<5 ~ e sin X Q to a very good approximation (even though e ~ 23°, the maximal difference is 
smaller than 0.3°, i.e. negligible in comparison with our experimental uncertainties). Making 
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the further simplification that the solar ecliptic longitude A Q increases linearly in time (i.e. 
assuming here that the eccentricity of the Earth orbit is e — 0), we have A Q = 2n(t — t )/T, 
with to the date of spring equinox. Combining this simple sinusoidal approximation for Sq (t) 
and Eq. (C.2) finally yields Eq. (2) of the main text. 



Appendix D. A simple model for the equation of time 

Following e.g. Ref. [1], a good approximation of the theoretical expression of the equation of 
time can be obtained in the following way. If the eccentricity of the Earth orbit were e = 0, 
and if the obliquity of the Earth were e = 0, one would have E — 0. We can thus expect that 
by calculating separately the small contributions E\ and E 2 of both the obliquity and of the 
eccentricity, and adding them, a good approximation of E is obtained: one basically expands 
E to the lowest orders in e and e. 

We first calculate the contribution E\ to E arising from the nonzero value of the 
obliquity. Here, we can neglect the ellipticity of the Earth orbit, and assume that the motion 
of the Sun around the Earth takes place on a circle, and thus, using Kepler's second law, at a 
constant angular velocity. The longitude A Q of the Sun along the ecliptic (see Fig. Aid) thus 
increases linearly in time as 

A = — (i-t ), (D.l) 

where T is the length of a year and to the date of the vernal (i.e. spring) equinox. However 
the mean Sun is a fictitious body that moves at constant angular velocity along the celestial 
equator, not the ecliptic. The difference between A Q and a , once converted to time via the 
correspondance 1° <H> 4 min, thus gives the obliquity contribution to the equation of time: 

Ek = (4 min/°) x (Aq - a Q ) (D.2) 

(for readability, from now on I shall drop the subscript 0). We thus have to express a as 
a function of A, and then use (D.l), to obtain E\(t). To find the relation between a and A, 
we introduce the equatorial frame {Oxyz) with origin at the center of the celestial sphere, 
the z axis pointing towards the celestial pole, and the x axis towards the vernal point V, and 
another, ecliptic frame (Ox'y'z 1 ) obtained from the former by a rotation of angle e around 
Ox. The coordinates of the Sun in the equatorial frame are (cos S cos a, cos S sin a, sin 6), 
and in the ecliptic frame (cos A, sin A, 0). Since the rotation matrix from the equatorial to the 
ecliptic frame reads 

10 \ 

cose sine , (D.3) 
— sin e cos e / 

one can relate (a, 5) to A. We obtain: 

tana = tan A cose, (D.4) 
sin S = sine sin A. (D.5) 

Therefore, using (D.4): 

A — a = A — arctan (cos e tan A) . (D.6) 

From this expression it may not obvious to see the time dependence (in particular the 
periodicity) of E\. We can obtain a better understanding (and a convenient expression) by 
noticing that cose is actually close to one (for e = 23.44°, we have cose ~ 0.9174). If 
one Taylor-expands f(x,A) = x — arctan(Atancc) around A = 1, one gets f(x, A) ~ 
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Figure Dl. Circular (gray) and elliptical (black) orbits around the Sun S with the same semi 
major axes a. The ellipticity of the ellipse is e = 0.3 for clarity. 



(1 — A) sin(2x)/2, and thus, putting everything together, we obtain our final expression for 

/ s d 1 — cose . f An . . \ _ „, 

El{t) ~2^ 2 sin (^(*-*o)J ■ (D.7) 

where d is the duration of a day. The dotted line on Fig. 4b shows Ex(t). 

We now turn to the calculation of E%, the contribution of the eccentricity of the Earth 
orbit: since its orbit is elliptic, the Earth does not move at constant angular velocity along its 
orbit, and the difference in angular position between the Earth and a fictitious body moving 
at constant speed with the same period gives the contribution E2 to the equation of time. 
Figure Dl shows the trajectory of the Earth E (black ellipse with the Sun S at one focus) and 
of the fictitious Earth E' having a circular orbit centered on S, with the same orbital period as 
E. From Kepler's third law T 2 /a 3 = Att 2 / (GMs un ), the radius r of the circular orbit is thus 
equal to the semi-major axis a of the elliptical orbit of E. Using the perihelia P, P' as the 
origins of angles, the polar angles defining the positions of P and P' are 8 and 8', respectively. 
We are interested in finding the difference in angular positions $ = & — 8, as a function of 
8, and then as a function of time. We will perform the calculation by keeping only first-order 
terms in the eccentricity e. 

The equation of the elliptical orbit of E reads [9] 

r E - fi^a (D.8) 
1 + e cos 8 

Moreover, Kepler's second law about the areal velocity implies that 

l -a 2 8> = *4r (D-9) 
2 T 

for the circular orbit, and 

\rtt> = %^ = (D ' 10) 



for the elliptical orbit, where we have used S^mpsc = nab with b = ayl — e 2 the semi-minor 
axis. 

Thus we get 

8' r 2 1 (1 - p 2 W 2 

t = ^-=L== ,\ L e) a , -l-2ecos0 (D.ll) 
8 a 2 Vl - e 2 (l + ecos6>) 2 

where the last approximation is valid to first order in e. Now, we have 
dz? d.0' 0' 

_ = _-!= l^-fecoB* (D.12) 



Measuring the eccentricity of the Earth orbit with a nail and a piece of plywood 



12 



and therefore 



i? = -2esin0. 



(D.13) 



Since we are keeping only first-order terms in e, we can replace 8 by 0' — 2ir{t — ti)/T in 
the above equation, where t\ is the time of perihelion passage. We finally get the following 
expression for E 2 : 



The dashed line on Fig. 4b shows E2(t). Combining (D.7) and (D.14), we find the expression 
(4) given in the text (solid line on Fig. 4b). 
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